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Abstract 

Due to the quantum correlation between the pair-produced D° and D° from the decay of the 
V>(3770), the time-integrated single and double tag decay rates depend on charm mixing amplitudes, 
doubly-Cabibbo-suppressed amplitudes, and the relative strong phase 5 between D° and D° decays 
to identical final states. Using 281 pb _1 collected with the CLEO-c detector on the ■0(3770) 
resonance, we measure the absolute branching fractions of D° decays to hadronic flavored states, 
CP eigenstates, and semileptonic final states to determine the relative strong phase, cos 5, of the 
K~ir + final state and to limit the mixing amplitude y. The results presented in this document are 
preliminary. 
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TABLE I: ST and DT yields for C = — 1 D°D° events, to leading order in x and y. 



When D° and D° mesons are pair-produced in e + e~ collisions with no accompanying 
particles (such as through the -?/>(3770) resonance), they are in a quantum-coherent C — — 1 
state. Because the initial state (the virtual photon]has J PC = 1 , there follows a set of 
selection rules for the decays of the D° and D° [2, IE IE HI IE IE 0- For example, both D° and 
D° cannot decay to CP eigenstates with the same eigenvalue. On the other hand, decay rates 
to CP eigenstates of opposite eigenvalue are enhanced by a factor of two. More generally, 
final states that can be reached by both D° and D° are subject to similar interference effects. 
As a result, the apparent D° branching fractions in this D°D° system differ from those of 
isolated D° mesons. Moreover, using time-independent rate measurements, it is possible to 
probe the D°-D° mixing parameters x = AM/T and y = Ar/2r, which are the mass and 
width differences between Dqp+ and Dcp-, as well as the relative strong phases between 
D° and D° decay amplitudes to any given final state. 

We implement the technique presented in Ref. |8(, where four types of final states are 
considered: flavored (labeled by / and /), CP+ eigenstates (S+), CP— eigenstates (£-), 
and semileptonic (£ + and t~). Event yields are functions of the number of D°D° pairs 
produced (denoted by Af, branching fractions (denoted by B) , the mixing parameters y and 
R M = (x 2 + y 2 )/2, and the amplitude ratio (f\D )/ (f\D°), whose magnitude and phase are 
denoted by rj and —Sf, respectively. We define Zf = 2 cos 5/ and give expressions for these 
yields in Table HI to leading order in x and y. 

Our analysis uses 281 pb _1 of e + e~ collisions, taken on the ^(3770) resonance, with 
y/s = 3773 MeV. The data were collected with the CLEO-c detector, which is a modification 
of CLEO III |E 0, 11, l3|, in which the silicon-strip vertex detector was replaced with 



a six-layer vertex drift chamber, whose wires are all at small stereo angles to the beam 
axis [13]. The hadronic final states we reconstruct are K~tt + (/), K + tc~ (/), K~K + (CP+), 
(CP+), K° s 7i°7i° (CP+), and K° s tc° {CP-). In the case of the two flavored final 
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states, K~ir + and K + tc~, both of these can be reached via Cabibbo-favored (CF) or doubly- 
Cabibbo-suppresssed (DCS) transitions. The strong phase between the CF and DCS d ecay 
amplitudes, 5k-k, is a source of ambiguity in some previous studies of D°-D° mixing [14]. 
We measure yields of both single tags (ST), which are single fully-reconstructed D° or 
D° candidates, and double tags (DT), which are events where both the D° and D° are 
reconstructed. We identify hadronic D candidates by their beam-constrained mass, M = 



^beam ~ Pc> and b y AE = E D - E hcam . 

We also measure semileptonic DT yields, where one D is fully reconstructed in one of the 
above hadronic modes and the other D is required to be semileptonic. We do not reconstruct 
semileptonic single tags because of the undetected neutrino. We also omit the DT modes 
where both D° and D° decay semileptonically. To maximize efficiency, we use inclusive, 
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partial reconstruction of the semileptonic D, demanding that only the electron be found. 
When the electron is accompanied by a flavor tag (K~tt + or K + tt~), we further require 
that the electron and kaon charges be the same, forming a Cabibbo-favored DT sample. 
Doing so reduces the dominant electron backgrounds, 7 — > e + e~ and tt° — > e + e~7, which 
are charge-symmetric. Such a requirement is unavailable for CP-eigenstate tags because 
they are unflavored. 

Efficiencies, backgrounds, and crossfeed among signal modes, are determined from Monte 
Carlo (MC) simulations. Following the least-squares procedure described in Ref. [l5j . we 
perform a fit to these efficiency-corrected yields to extract the free parameters listed above. 
We assume that Kg is a purely CP-even eigenstate and that CP violation in D° decays is 
negligible. In Table HTl we show the preliminary results of the data fit. Because the precision 
of the world average for r 2 KlT far exceeds our determination [3, 17 , [l8| , we constrain this 
parameter to be (3.74 ± 0.18) x 10~ 3 in the fit. The \ 2 is 15.7 for 20 degrees of freedom, 
and only statistical uncertainties have been included. Systematic uncertainties are being 
evaluated, and it is expected that they will be of similar size. The value of £>(P° — > K^tt ) 
shown in Table ITT1 is equivalent to and correlated with the so-called "single tag" measurement 
in Ref. 3 of (1.212 ± 0.016 ± 0.039)%, which is based on the same dataset as the current 
analysis but makes use of independently-performed measurements of y 16|| and M jib] ], 
whereas we allow both of these parameters to be determined by our fit. 

As discussed in Ref. systematic effects that are correlated by final state, such as 
mismodeling of tracking or 7r° reconstruction efficiency, cancel in the DCS and mixing pa- 
rameters. However, one important source of uncertainty is the quantum-number purity of 
the reconstructed CP eigenstates. Peaking backgrounds to CP eigenstates may come from 
flavored decays or CP eigenstates of the opposite eigenvalue. Therefore, the size of the 
simulated background, which assumes uncorrelated decay, may differ from reality because 
the quantum correlation modifies the rates of each of these processes in a different way, and 
a systematic uncertainty can be assigned based on the fit results. 

Also, the purity of the C — — 1 initial state may be diluted by radiated photons, which 
would reverse the C eigenvalue. We limit this effect by searching for DT modes with same- 
sign CP eigenstates (such as K~K + vs. 7r + 7r~). These decays are forbidden for C — — 1 
but are maximally enhanced for C = +1. Including these yield measurements (all of which 
are consistent with zero) and fitting all the other yields to a sum of C = — 1 and C — +1 
contributions, we find no evidence for C = +1 contamination — the C = +1 fraction of the 
sample is 0.06 ±0.05 (stat.) — and we observe no significant shifts in the fitted parameters. 

In summary, using 281 pb _1 of e + e _ collisions produced on the -?/>(3770) at CLEO-c, 
we have searched for D°-D° mixing and made a first measurement of the strong phase, 
5ktt- We expect future improvements with the addition of more CP eigenstate modes, more 
■^(3770) data, and higher-energy data with P°P°7 events, where the D°D° pair is a C = +1 
eigenstate. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. D. Cronin-Hennessy and A. Ryd thank the A. P. Sloan 
Foundation. This work was supported by the National Science Foundation, the U.S. De- 
partment of Energy, and the Natural Sciences and Engineering Research Council of Canada. 
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Parameter Fitted Value PDG [16] 
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TABLE II: Preliminary results from the data fit, with r\ v constrained to be (3.74 ± 0.18) x 10" 
Uncertainties on the fit results are statistical only. 
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